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1. Introduction 



The equation of motion for the complex sl{n)^^^ Toda field theory can be written 



2 

m 



0(f) = — — ctj exp ijSaj ■ (j). (1-1) 

The field 4>{x, t) is an n — 1-dimensional vector: an element of the Cartan subalgebra of the 
finite Lie algebra sl{n). The inner products are taken with respect to the Killing form of 
sl{n) restricted to the Cartan subalgebra. The ctj's, for j = 1, 2, . . . , n — 1 are the simple 
roots of sl{n); ao is the extended root (minus the highest root of sl{n)). The fact that the 
extended root is included in the sum, distinguishes the affine theories from the non-affine 
ones. The Q!-,'s are linearly dependent: 

n 

E«i = 0- (1-2) 

i=i 

The constant m sets an arbitrary mass scale and /? is a coupling constant. In the complex 
theories P is real, whereas, in the real theories /? is purely imaginary. 

The real theories, when ^5 = z/? e R, are well understood, both in the classical and 
quantum regime (see for example [1]). The spectrum consists of n — 1 particles of mass 

Tra 

ma = 2m sin — , a = 1, 2 . . . , n — 1. (1-3) 

n 

In the quantum theory the spectrum is preserved, except for an overall renormalization of 
the mass scale m. It is known from a Feynman diagram calculation to one- loop [1], that 
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An n 



(1.4) 



(In the following we shall use "hats" to denote exact quantum masses.) The complex 
theories, when /3 G M, have a much more complicated spectrum, since they admits kink or 
soliton solutions. This is well-known when n = 2, for which the real case is the sinh-Gordon 
theory and the complex case is the sine-Gordon theory. Soliton solutions were first written 
down in ref. [2] for the affine sl{n) and d4 theories; recently, soliton solutions have been 
found for all the affine Toda theories [3] . 
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For the sl{n) theories the general A'"-sohton solution is constructed in the following 
way. To each soliton one associates the data {cr, A, a, where a and A are real parameters 
satisfying 

:r((7, A, a) = 0, (1.5) 



where the characteristic polynomial is 



Tva 



T{a, A, a) = a-\- 4m^ sin^ — 



n 



(1.6) 



a is an integer in the set {1,2, ...,ri — 1} and ^ is, for the moment, an arbitrary complex 
parameter. To each soliton, say the p*"^, one associates the n functions 

n 



^^f'ix, t) = apX - Xpt + ^apj + ^p, j = 1,2,..., n, (1.7) 
and to each soliton pair, say the p^'^ and q^'^, one associates the "interaction function" 

(1.8) 



^{cFp + C^g, Ap + \q, Up + ttg) 

The general A'^-soliton solution to the equations of motion is 
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Mi=0 Hn=0 

So the one-soliton solution is 



p=i ^<p<q<N 



(1.9) 



1 



i=i 



1 + exp I crx — At + ^ H j 



n 



with 



cr^ - A^ = 4m^sin^ — . 

n 



(1.10) 



(1.11) 



The soliton is a kink, whose centre-of-mass is at cr~^[At — Re^], moving with velocity A/u 
and having characteristic size . The parameters a and Im^ determine the topological 
charge of the soliton. The topological charge is defined to be 



2n 



d(f) 



t = — / dx — — , 



dx 



(1.12) 



from which one readily verifies that (1.10) has a topological charge which a weight of the 
a*^ fundamental representation of sl{n), where the exact weight is determined by Im^. 



The classical masses of the solitons are easily determined by explicit calculation from 
the Hamiltonian: 



One finds 

Mo = -T^sm— , a= l,2,...,n- 1. 1.14) 

Notice that these masses are real, despite the complex form for the Hamiltonian. More 
particularly, notice that they are proportional to the fundamental masses of (1.3): Ma oc 
nia- Both these facts were explained in [2]. 

In ref. [4] the quantum theory of the solitons was investigated. A form for the soliton- 
soliton S'-matrix was proposed, which agrees with the WKB quantization of the classical 
scattering theory, discussed above, in the semi-classical limit. In general, it was shown 
that the S'-matrix represents a non-unitary quantum field theory, but that there was a 
unitary region depending on the coupling constant, for which, in addition, the bootstrap 
closes on a finite set of states corresponding to solitons transforming in the fundamental 
representations of sl{n). Central to the proposal in ref. [4] for the quantum theory of 
the solitons, is the idea that the ratios of the soliton masses should not be altered in the 
quantum theory: Ma oc nia- It is clearly of interest to calculate the quantum corrections 
to the soliton masses directly, in order to verify the truth of this fact. 

There exists a standard method for computing the first quantum corrections to soliton 
masses in a 1 + 1-dimensional field theory: one basically sums the zero-point energies of 
the fiuctuations around the soliton solution. Obviously, such a sum is divergent, but these 
divergences may be removed, in the standard way, by renormalization. Moreover, there is 
a piece of folklore for 1 -|- 1-integrable theories that says that the masses of solitons obtained 
in this way are exact and there are no higher-order corrections. This is appears to be true 
for the sine-Gordon theory (the case for g = sl{2)), where the exact soliton mass is 

8m 2m ,^ ^ ^, 

M=^--, (1.15) 

in our conventions — the first term being the classical contribution from (1.14), and the 
second term being the first quantum correction. This piece of folklore can, to some extent, 
be justified by noting that at some level the theory can be cast, via a transformation to 
action-angle variables, into a set of harmonic oscillators — albeit of infinite number — and 
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the saddle-point calculation of the soliton masses should, therefore, be exact. Although 
this argument has a certain verisimilitude, it is not clear how general the principle is; in 
particular, does it apply to Toda theories, beyond the sine-Gordon theory? 

As a by-product of our calculation of the quantum corrections to the soliton masses, we 
are able to prove that the soliton solutions are classically stable; an issue that is important 
given the fact that the Hamiltonian (1.13) is complex. 

2. Classical Stability 

In this section, we consider the question of the classical stability of the one-soliton 
solutions. This question is closely related to the calculation of the quantum corrections to 
the soliton masses, as we shall see. 

The standard way to approach the question of stability, is to consider the effect of 
adding a small perturbation 77 to the one-soliton solution (1.10), which we denote 0. It 
is convenient to take the soliton to be stationary with the centre-of-mass at the origin: 
A = 0, cr = nia and Re^ = 0. To first order in 77 the equation of motion (1.1) becomes 

n 

Drj + '^ji'^j ■ v) 6xp i(3aj ■ (f) = 0. (2-1) 

Now consider a perturbation with a time-dependence of the form r]{x, t) = fj{x) exp iut, so 
that 

Vfj = u^f], (2.2) 
where we have defined the following differential operator: 

q2 ^ 

V — + '^j ® 6xp i(3aj ■ 4>{x). (2.3) 

The question of stability now boils down to showing that the spectrum of P — for bounded 
eigenfunctions — is real and positive; hence, the frequencies v are real and small perturba- 
tions to do not diverge. Our strategy for proving this will be to find the exact spectrum 
of V. 

First of all, the spectrum of V is real, since it is actually a hermitian operator, with 
respect to some inner product. To show this, we first notice that the one-soliton solutions 
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— for some particular choice of Im ^ — actually satisfy a reality condition of the form 



4)*{x) = -M4){x), (2.4) 
where M acts as a Z2 symmetry of the roots aj: 

Maj=am(j), rn'^{.j)=j- (2.5) 
This is explained in detail in refs. [2,5]. From this it follows that 

= MVM. (2.6) 
Hence, V is hermitian with respect to the inner product 

/oo 
dxp-Mg, (2.7) 
-00 

for functions f{x) and g{x) in the Cartan subalgebra of sl{n). 

We now turn to the determination of the spectrum of bounded cigcnfunctions of V. 
The eigenvalue equation (2.2) can be viewed as a multi-component Schrodinger problem 
for the potential 

n 

The spectrum has contributions from two sources: the bound states and the scattering 
solutions. The former have a discrete spectrum and the latter a continuous spectrum. 
There are also eigenfunctions which are not bounded, but these have no relevance for the 
discussion of stability and for the quantum mass corrections. 

In general the scattering solutions have the following asymptotic behaviour 

lim f]{x) — Aexp ikx + B{k) exp —ikx 

(2.9) 

lim fi{x) = C{k) exp ikx, 

X—KX) 

where the incoming, reflection and transmission coefficients. A, B{k) and C{k), are Cartan 
subalgebra- valued. Remarkably, the potential in (2.8) is "reflection- less" , so that B{k) = 0. 
This can be shown directly using the following argument. In the introduction, we wrote 
down an expression for the A^'-soliton solution. Consider a two-soliton solution where the 
first soliton is 0, the static one-soliton solution that we started with, and the second soliton 
— the "probe" — will act as a small perturbation of the first. To this end, we take a2 = ik, 
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02 = b, X2 = —iv and treat exp ^2 as a small parameter; in which case to first order in 
exp ^2 we deduce 

1 + uj^^ exp(maa; + ^ + 7) 



fib{k\x) = ^oij 



uj'^^ exp ikx, 



(2.10) 



1 + ijo°-i exp(maa; + ^) 

where we have introduced the primitive n^^ root of unity. In the above, 7 is the 
interaction parameter (1.8) which gives in this case 

— 2imak 



exp 7(/c) 



(2.11) 



The expression in (2.10) is the exact scattering solution to the Schrodinger problem (2.2) 
with eigenvalue = + ml and 

n n 

B{k) = 0, A= ^u^^aj, C{k) = exp 7(fc) J^cu^^a^. (2.12) 

Moreover, the set of solutions fjb{k; x), for 6 = 1, 2, . . . , n — 1 and /c e M, forms a complete 
set of scattering solutions, since the vectors Yl^=i ^^^OLji for 6 = 1, 2, . . . , n — 1, span the 
Cartan subalgebra. Hence, (2.8) is a reflection-less potential. 

The beauty of dealing with a reflection-less potential is that it allows for a simple 
determination of the bound-states as well. These solutions occur for values of k for which 
the transmission coefficient has a zero, and which are, in addition, normalizable. It is easy 
to see, from the explicit expression for exp7(/c) in (2.11), that C{k) has a zero when 

7r(a — h) 



ml + ml- ml_^ 

k = = —imi, cos 

zirrin 



n 



(2.13) 



The exact solution for the bound-states follows from (2.10) and the particular value for k 
in (2.13): 



Vb 



co^^ exp 



a-, 



rribX cos 



77(0—6) 



1 + u;"-? exp(maX + ^) ' 



which has a frequency 



. 7r(a — b) 
V = rrib sm 



n 



(2.14) 



(2.15) 



The asymptotic limits of the above solution are 



lim r]{x) = ajUJ^^ exp 



X—I- — 00 



„ . nb 7r(a — b) 

2mx sm — cos 

n n 



lim ri{x) = 'S^ a-jUJ^'^ "^-^ exp 

;r — ►no f ^ 



. 7r(a — b) Tzb . 

2mx sm cos ^ 

n n 



(2.16) 



Of course, not all these solutions have the right asymptotics to be bona-fide bound-states of 
the Schrodinger problem. The bone-fide bound-states must be normalizable which means 



lim ?7(a:) ~ exp |a;, (2-17) 

for some real non-zero constants k,^ ^ 0. The values of b which lead to bona-fide bound- 
states are 

71/ fX Tl 

1 <b < a or — < b < — + a, fora< — 

n , 72 , / 1 r ^ (2.18) 

a < < — or a < < n — 1, tor a > —. 

2 2 - - 2 

The solution with b = a has not been included since this corresponds to the zero-mode of 

V given by 

Vaix) = > : ; -r. 2.19 

J^i 1 + exp (m^a; + 

The appearance of this zero-mode was only to be expected since it is proportional to 
d4>/dx, and is due to the freedom to shift the centre-of-mass of the soliton. 

In addition to these bound-states there are bounded solutions which must also be con- 
sidered. Such eigenfunctions have a constant asymptote, and hence are not normalizable. 
They appear only when n is even and are given by (2.14) when b = ^or6 = a-|-^, if 
a < |, and 6 = a - | or 6 = |, if a > |. 

From the preceeding analysis, we conclude that the spectrum of bounded eigenfunc- 
tions of T> is real and positive, except for the zero-mode which reflects the freedom to move 
the centre-of-mass of the soliton. Hence the one-soliton solutions are classically stable to 
small perturbations. 



3. Quantum Corrections to the Masses 



In this section we calculate the first quantum corrections to the soliton masses. We 
employ the semi-classical WKB method described in [6] for the kinks, or solitons, of the 0^ 
and sine-Gordon theories — see also the review in [7] . The method is completely standard; 
however, we feel that the Toda theories are sufficiently more complicated than these other 
examples that it is worthwhile including all the calculational details. In particular, the 
question of boundary conditions is not completely standard. 



7 



The dimensionless expansion parameter in our theory is jvn? ^ and so the expansion 
in h coincides with the weak-couphng expansion. From now on we set ^ = 1. The WKB 
method of [6] is straightforward to apply in the present situation since the one-sohton 
solution is, in its rest- frame, time-independent. The idea is to compute the zero-point 
energy of the small oscillations around the classical solution. The sum over all the modes 
can be done when an infra-red regulator is introduced; this is most easily achieved by 
putting the theory in a box with rigid boundary conditions. The zero-point energy of the 
vacuum must then be subtracted. The resulting expression is independent of the size of 
the box, L, as L — ^ oo; however, the final expression has an ultra-violet divergence. This 
divergence is not a problem of the soliton solution per se, rather, it is just a manifestation 
of the fact that the bare-mass needs to be renormalized. This is achieved simply by normal- 
ordering the Hamiltonian. The extra correction removes the divergence and the resulting 
finite residue which remains is then the mass correction. 

The small fluctuations which contribute to the mass shift come from either the bound- 
states or the scattering solutions of the linearized equation of motion around the soliton 
solution: the Schrodinger problem (2.2). 



3.1 Contributions to the Mass from the Bound-States 



Each of the bound-states fjb{x), with b in the range (2.18), contributes an amount 
to the mass shift, where u is given by (2.15). In addition to this, the bounded solutions 
— which only occur for n even — also contribute to the mass shift. However, we shall 
not include the states (2.14) with b = a + fora< ^, and 6 = a— ^, fora> ^, which 
correspond to the k = solutions of (2.10), and whose contribution to the mass shift will 
be included in the contribution from the scattering states discussed in the next subsection. 
So the contributing modes of the form (2.14) are 



1 < b < a or — < b < — ha, fora< — 
2 ~ 2 2 



a 



2 ~ 2 



l<b< 



n 



or 



or 

n 



a < b < n — 1, 
<b <n-l, 



for a > 



for a = |. 



n 



(3.1) 



The contribution to the mass shift from these states, for a < ^, is 



AMa{l) = m 



. Trb . 7r(a — b) 

> sm — sm 

^-^ n n 

6=1 



< 2+« 



E. irb . 7r(a — b) 
sm — sm 
n n 



(3.2) 



(Notice, that the expression is also vahd for a = ^, since in that case the term with = f 
does not contribute.) It is straightforward to evaluate the sums involved, giving 

It is not difficult to show that the cases when a > ^ give the same results as (3.3), as a 
function of a. 

3.2 Contributions to the Mass from the Scattering States 

In this subsection, we calculate the contribution of the scattering solutions of (2.2) 
to the mass correction. The idea is to sum the zero-point energies of all the scattering 
modes. Such a sum will, of course, be infra-red divergent, so first of all it is necessary to 
introduce some finite boundary conditions so that the modes become discrete and therefore 
enumerable. Usually, periodic boundary conditions are chosen; however, they are not 
appropriate here because j{k) in (2.11) is not purely imaginaginary. It turns out that the 
appropriate boundary conditions are •fi{x = 0) = 77(0; = L) = 0, where the centre-of-mass 
of the soliton lies somewhere inside the box, and L is much larger than the size of the 
soliton: L ^ ^. The solutions satisfying these boundary conditions are 

Vbikp-, x) - fjb{-kp; x), (3.4) 

where fjb{k; x) is the solution in (2.10), with 

kpL + pb{kp) = 7rp, kp>0, p e Z. (3.5) 

where Pb{k) — Im7(A;) and 7(/c) is defined in (2.11). 

The quantum correction to the soliton mass is then obtained by taking the zero-point 
energy of the modes around the soliton solution and subtracting the zero-point energy of 
modes around the vacuum, which satisfy (3.5) with p = 0. That is 

n— 1 ^ V 

AM,(2) = i E E W^P + ™6 - \/[kp + Pb{kp)/Lf + , (3.6) 
b=l kp>0 ^ ^ 

where the sum over kp is the sum over distinct solutions of (3.5). Since L will eventually 
be taken to infinity, we can expand in L~^, the leading term in the sum being 

-L-^^l^ = -L-^,,(k)'j^, (3.7) 
^/k^ + ml dk 

9 



with eb{k) = •\/^^~+m|^. Now, we take L — > oo in which case we can replace the sum over 
kp by an integral over k: 



p>0 



dk 
2^ 



+ 0{L-'] 



(3. 



where we have used the fact that the integrand is an even function of k. Therefore, the 
mass correction is 



6=1 -^-^ " 



deb{k) 
dk 



^ n-1 r ^oo 

-JzJ2 Pb(^'>^b{k)\'^^+ dkeb 
6=1 L 



{k) 



dpbjk) 
dk 



(3.9) 



Now as k ^ ±oo, ei,{k) \k\ and 



{2ml + - m^_(, - m^+(,) + ^■ 



(3.10) 



and thus the first contribution from (3.9) is 



^ n— 1 ^ n — 1 



m 



a—h 



(3.11) 



6=1 6=1 
The sum is straightforward to perform, the result being a contribution 



27r 



(3.12) 



to the mass. As it stands, the second contribution in (3.9) is not well defined because the 
integral has a logarithmic divergence. The problem is not due to any special nature of 
the soliton solution, indeed it is to be expected, being due to the fact that we have not 
renormalized the theory. We now pause to consider this aspect in more detail. 

Mercifully, the renormalization process in a two-dimensional field theory is straight- 
forward. The divergences can simply be removed by normal-ordering the Hamiltonian. 
Working to lowest order in 0^ and introducing an ultra-violet cut-off A 



: exp i(5aj • := exp i(5aj 



n-1 



0=1 



where 



^ _ dk 1 



m 



2' 



(3.13) 



(3.14) 
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and Ca is the a^^ eigenvector of Y^j^i'^j ® eigenvalue (ma/m)'^. Hence, to lowest 

order in the renormalized Hamiltonian is equal to 

with 

n-l 

= i(m/3)2 J^iCa ■ «,)'A«. (3.16) 
0=1 

The presence of the additional terms clearly changes the energy of the soliton solution. 
These changes must be added to the quantum correction to the soliton mass. It might 
be thought that we should now consider the one-soliton solution for the renormalized 
equations of motion; fortunately, this is unnecessary since the correction due to change in 
the form of the solution does not occur at the lowest order in (3'^, precisely because the 
soliton solution satisfies the "bare" equations of motion (1.1). 

Taking into account the counter-term from the renormalization, the final expression 
for the mass correction for the a^^ soliton from the continuum solutions is 



, man v-^ f dk 



dk 



1 

— —2 / dx [ex.p{i(3aj • 0) — l] dm 

^ 3 = 1 



(3.17) 



—00 



We now pause to calculate the integral 

/oo 
dx [exp{iPaj • ^) — l] . (3.18) 
-00 

Using the explicit form for the stationary one-soliton solution, one finds 

exp(zpQ!,- • 0) — 1 = — — 7T. (3.19) 

' Vm7 [1 + u'^J exp{max + 0] 

The integral (3.18) may now be calculated explicitly yielding 

"rngV- 
. m ) 



V m 7 mo [1 + (jj"-^ exp{maX + ^)] 

11 



nig 



(3.20) 



The important point about the resulting expression for the integral is that it is in- 
dependent of j. The contribution to the mass correction (3.17) from the renormalization 
counter-term can now be simplified: 



-j^ POO 

— — ^ / [exp{iPaj • 0) — l] dm 



n— 1 n 
6=1 j^l 

fmi,\'^ dk 



(3.21) 



n— 1 „ „A 

--SOL 



A 47r ypqPm^ 



where we used the fact that (i, is an eigenvector of X^^^i ocj aj of eigenvalue {rrib/m)" 
The expression for the mass correction is now 

n-l 



ml 



V + ™^ (m2 + - ml^,Y + 4m2/c 



(3.22) 



Using the fact that 



n-l 



6=1 



and 



n-l 



6=1 



V m 



Tr 



/ n 



(3.23) 



(3.24) 



one can easily verify that for large the integrand behaves as k~'^ and so, as required, 
the logarithmic divergence is cancelled by the renormalization counter-term. We can now 
let the ultra-violet cut-off A tend to infinity. 

The remaining contribution is given in terms of a convergent integral: 



n 



AMa{2) = ma [-— +/(a,n) 



(3.25) 



where 



mf 



{ml + ml- ml^^Y + Amlk"' 



1 



(3.26) 
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Unfortunately, we have not managed to evaluate the integral /(a, n) analytically; except 
when n = 2 in which case the integrand is zero. However, by evaluating the integral 
numerically for a range of a and n we find, to a very high degree of accuracy, that, firstly, 
the integral is independent of a, and secondly it is equal to the following functional form 

[8]: 

/(a,n) = |-!T.'^ '"^^'^'^ (3-27) 

^ ' ^ I -^cot- n e Even. ^ ^ 

We can put the two contributions from the bound-states and the continuum together 
to arrive at our main result: 

(Tl 1 TT \ 

-— + -cot- . (3.28) 
/TT / n J 

Notice that the final result is valid for n even or odd: the asymmetry has cancelled out 
on adding the two contributions. 



4. Discussion 



We have succeeded in calculating the first quantum corrections to the solitons masses 
in complex sl{n) affine Toda field theories; and as a by-product proving that the solitons are 
classically stable. Although the calculations were somewhat lengthy the result is simple: 



Ma = 2nma 



1 



1 



1 



An An n 



(4.1) 



where the first term is the classical mass. For the sine-Gordon theory it is thought that 
the result is exact and there are no higher order corrections; this example is recovered by 
setting n = 2, giving 

M - ^. (4.2) 



where P' is the characteristic "renormalized" coupling 

l-/32/47r' 



(4.3) 



In the sine-Gordon theory it is well-known that the fundamental particle is, in the quantum 
theory, a soliton anti-soliton bound-state. However, it has been conjectured in ref. [4] that 
the same is true in the more general theories. The soliton-soliton ^'-matrix written down 
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in ref. [4] has poles corresponding to soliton bound states; in particular there are n — 1 
distinct sets of scalar states, carrying zero topological charge. The ground states of these 
sets, have mass 

ma = 2MaSm^, a = 1, 2, . . . , n - 1, (4.4) 

where Ma is the soliton mass and A is a coupling constant which is known to be related 
to P as [4] 



A = ^ + A, (4.5) 



471 

J' 

where A is 0(1), but was not determined in ref. [4]. By using (4.1) and (4.5) we can 
expand (4.4) to first order: 



m„, = rur. 



' 4n n Att Att ' ^ ' 



(4.6) 



But (4.6) should be compared with the one-loop formula for the masses of the fundamental 
particles in (1.4) (with $ = if3). They are consistent if A = —1 so 

^=^' (4-7) 

to this order in /3, where (5' is the shifted coupling appearing in the sine-Gordon theory 
(4.3). 

It would be interesting to investigate the question as to whether the masses in (4.1) 
are actually exact. On way to do that would be to repeat the analysis of ref. [6], and 
calculate the corrections to the masses of the breathers. 

I would like to thank M.R. Niedermaier for finding the analytic expressions for the 
integrals I{a,n). I also have pleasure in thanking J.M. Evans for contributing to the 
investigation of the classical stability of the solitons in Section 2. Finally, I acknowledge 
the support of Merton College, Oxford, through a Junior Research Fellowship. 
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